Quasi-Local Linear Momentum in Black-Hole Binaries 
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We propose a quasi-local formula for the linear momentum of black-hole horizons inspired by the 
formalism of quasi-local horizons. We test this formula using two complementary configurations: 
(i) by calculating the large orbital linear momentum of the two black holes in an unequal-mass, 
zero-spin, quasi-circular binary and (ii) by calculating the very small recoil momentum imparted to 
the remnant of the head-on collision of an equal-mass, anti-aligned-spin binary. We obtain results 
consistent with the horizon trajectory in the orbiting case, and consistent with the net radiated 
linear momentum for the much smaller head-on recoil velocity. 
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I. INTRODUCTION 



The major breakthroughs in numerical relativity that 
occurred two years ago [if, |3, Q not only allowed numer- 
ical relativists to faithfully simulate the inspiral, merger, 
and ringdown of arbitrary black-hole-binary configura- 
tions, but also inspired new developments in mathemat- 
ical and numerical relativity and astrophysics. These in- 
clude exploring cosmic censorship , modeling the hori- 
zon spin direction in non-axisymmetric spacetimes 
modeling black- hole — neutron-star binaries @ , modeling 
the collapse of supermassive stars Q, and finding recoil 
velocities in astronomical observations [i,ai3. These 
breakthroughs have also had significant consequences for 
the data analysis of LIGO and other ground based ob- 
servatories that are attempting to directly detect gravi- 
tational waves [ll|, [H, [ll] . 

The recoil velocities acquired by the remnant of the 
merger of black-hole binaries has many interesting as- 
trophysical consequences [l4| . particularly since spin- 
ning black holes can accelerate the merged hole up to 
4000km [IB], large enough to eject the remnant from 
the host galaxy. As was noted in [3, [H, [itI , the spin 
contributions to the recoil velocity are generally larger 
than those due to the unequal masses, and, in particular, 
the spin component in the orbital plane has the largest 
effect [ll] , leading to a maximum recoil velocity of about 
4000km The first study of generic black-hole- 

binary configurations (i.e. a binaries with unequal com- 
ponent masses and spins, and spins not aligned with each 
other or the orbital angular momentum) was described 
in Ref. [3|, and, based on the results of that study, a 
semi-empirical formula to estimate the recoil velocities 
of the remnant black holes was proposed, finding recent 
confirmation in [l^, [H, [l^ . 

In all the above papers the computation of recoil ve- 
locities was made by directly extracting the radiated lin- 
ear momentum of the system at large distances from the 
remnant black hole (asymptotically approaching future 
null infinity J^"*") [20|. Here we are interested in an al- 



ternative measure of linear momentum; one in terms of 
quantities defined local to the horizon. Such a formula 
would provide an independent computation of recoil ve- 
locities (and thus provides a further consistency check 
for numerical simulations) and would allow for the com- 
putation of the linear momentum of the individual holes 
while they orbit each other before the final plunge and 
merger. The starting point of our analysis is the quasi- 
local description of black hole horizons as provided by the 
isolated [Uj, dynamic al 22] . and trapping horizon [2^ 
formalisms; see [2J, [23, |2a| for reviews. These techniques 
have been successfully applied to black-hole binaries to 
extract the mass and spins (magnitude and direction) 
of the individual holes and the merger remnant, as well 
as measure the precession of the spin direction and spin 
flips [^ [27I [28l . [29| . In this paper we propose, and nu- 
merically test, a formula inspired by the theory of quasi- 
local horizons to compute the linear momentum of a black 
hole. 



II. QUASI-LOCAL LINEAR MOMENTUM 

In standard classical mechanics and field theory, con- 
served quantities are defined, using Noether's theorem, as 
the generators of symmetries of the action. Thus, angular 
momentum is defined as the generator of rotational sym- 
metry, energy is the generator of time translations, and 
linear momentum the generator of spatial translations. 
This approach has been successfully used to calculate the 
energy and angular momentum of isolated and dynami- 
cal horizons [SO, [Mi [13, [1^ ■ For classical mechanics and 
field theory in flat Minkowski space, there is a natural 
way to define translations and rotations using the sym- 
metries of the flat background metric. The situation is 
more complicated for general relativity on a curved space- 
time manifold when we typically do not have the luxury 
of being able to use a reference background metric; this is 
precisely the situation for numerical relativity black-hole 
simulations. To even speak about conserved quantities 
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in this context, one therefore needs to assume the exis- 
tence of an (at least approximate) rotational symmetry 
for angular momentum, and similarly a preferred time 
translation vector field at the horizon for energy. Since 
the spacetime in the vicinity of a black-hole binary is not 
translationally invariant, it is clear that the Hamiltonian 
calculations cannot be simply extended to define linear 
momentum for the individual black holes. Because of this 
conceptual difhculty, we follow a different approach for 
linear momentum: we shall take the Hamiltonian calcula- 
tions for angular momentum and energy only as motiva- 
tions for our definition of linear momentum. As we shall 
see, while our results are very promising, there are open 
issues remaining; this work should therefore be seen as a 
preliminary investigation to be followed up with further 
mathematical and numerical work. 

Motivated by the analysis of angular momentum in 
[3^ . let us start with the momentum constraint on a 
spatial slice E of the spacetime with extrinsic curvature 
Kab and 3-mctric "fab- 

D''Pab=0 where Pab = Kab - K-fab ■ (1) 

Here Da is the derivative operator on S compatible with 
7ab. Let be a vector field tangent to E which looks 
asymptotically like a unit translation, and jab is the 3- 
metric. Let us first consider the case when there is a 
single apparent horizon 5, considered to be the inner 
boundary of S, and let 6*00 be the sphere at spatial in- 
finity. We are interested in the linear momentum along 
Contracting the momentum constraint with and 
integrating by parts, we get 

(f-(f PabCd^S" f P-'Caabd^'V, (2) 

where is the Lie derivative with respect to We 
have assumed asymptotic flatness at spatial infinity with 
Kab = 0{r~^) and jab = Sab + 0{r~^) at large distances 
from the source. Since our focus here is near the horizon, 
we shall not spell out any further details of the asymp- 
totic behavior. 

We would like to interpret Eq. ([2]) as a balance law for 
linear momentum. It should be noted that the analog of 
this equation for a rotation = ip'^ leads to the standard 
expressions for angular momentum at infinity and at the 
horizon [U ■ As expected, the surface integral at infinity 
is the usual ADM momentum associated with the entire 
Cauchy surface, and we would like to identify the surface 
term at S with the linear momentum of the black hole: 

Pf^ ^^j^{Kab-Kjab)Cd^SK (3) 

The factor of Stt is chosen to recover the ADM linear mo- 
mentum at 5*00 (and we have chosen units with G — 1). 
The right had side of Eq. ([2]) is interpreted (up to a fac- 
tor of Stt) as the "flux" of linear momentum between the 
horizon and infinity; it would vanish if is a Killing vec- 
tor of the 3-metric. We do not expect this flux term to 



vanish generically. However, it does vanish for maximal 
slices {K — 0) if f is a conformal Killing vector of the 
3-metric jab- In particular, the flux vanishes for the con- 
formally flat Bowen-York data for a single boosted black 
hole if we take to be any of the coordinate basis vec- 
tors d/dx^d/dy^d/dz. Thus, in this case agrees with 
the ADM linear momentum. For Bowen-York data with 
N multiple boosted spinning black holes, if we denote the 
i*'' apparent horizon as S**, the same argument as above 

shows that P^^ •* is additive: 

N 

pADAI^Y^PP- (4) 

Thus we see that, for Bowen-York initial data and, some- 
what more generally, for maximal conformally flat spatial 
slices, Eq. ([3]) can be considered as a satisfactory formula 
for quasi-local black hole linear momentum. Moreover, 
the derivation of the balance law Eq. ^ did not any- 
where use the fact that 5" is a marginally trapped surface; 

the balance law is valid for any inner boundary. In par- 

(S) 

ticular then for maximal conformally flat slices, P^ can 
be used to measure the linear momentum contained in 
closed regions in E. 

In the general case during an evolution, we will have 

neither conformal flatness nor maximal slices. The lin- 

(s) 

ear momentum flux will not vanish, and P^ will not 
agree with p^^^ , This fact, by itself is not necessarily 
a problem. The most obvious problem with interpreting 
as linear momentum then is the choice of at the 
horizon; what should it's direction be, and how should 
it be normalized? In the regime when the black holes 
are sufficiently far apart and the orbit is varying slowly, 
a possible choice for is (spatial projection of) the he- 
lical vector tangent to the orbit. Alternatively, the cor- 
rect method might be to find an approximate translation 
Killing vector field on E in the vicinity of the horizon by 
adapting the Killing vector finding methods described 
in 27] or 28]. In both cases it is not clear what the 
correct normalization of should be. If E is the quasi- 
local horizon energy corresponding to the particular time 
translation vector field used [sl], then it is not clear if 
there is a sense in which {E,P^) can be viewed as a bona 
fide energy-momentum four vector; for example, can we 
prove E'^ = P^ + M"^ where M is the horizon mass? Some 
of these issues could be studied using approximate initial 
data of the kind constructed in [s^ . 

Further study of these issues will be left to future work 
and in this paper, we will work with being one of the 
coordinate basis vectors in the coordinate system used in 
the numerical simulation; this will yield the three com- 
ponents of linear momentum. This is obviously gauge 
dependent, but we shall show that there exist suitable 
gauge choices in which the linear momentum passes ba- 
sic but non-trivial consistency checks. The situation is 
similar to what was observed for quasi-local horizon an- 
gular momentum in [sj . The true gauge invariant angular 
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momentum requires an accurate calculation of the axial- 
symmetry vector field on the horizon. However, in 
some cases it might be possible to calculate the com- 
ponents of angular momentum by using the rotational 
vector fields constructed from the coordinates used in 
the simulation; but this is by no means guaranteed. The 
same situation presumably holds for linear momentum as 
well. 



III. NUMERICAL TECHNIQUES AND 
RESULTS 

We use the puncture approach [s^] along with the 
TwoPuNCTURES [111 thorn to compute initial data. 

We confirmed that Eq. Q yields an accurate evalu- 
ation of the horizon momentum on the initial slice for 
Bowen-York data. In our tests we found agreement be- 
tween the momentum parameters and the calculated mo- 
mentum to better than 1 part in 10^ both for single 
boosted black holes with momenta in the range < 
P/M < 10 and for orbiting black-hole-binary datasets 
with coordinate separations as small as 0.25M and spins 
as large at a/m = 0.84. 

We evolve these black-hole-binary datasets using the 
LazEv '39*1 implementation of the 'moving puncture ap- 
proach' which was independently proposed in 0|, Q. In 
our version of the moving puncture approach we re- 
place the BSSN [10, Hi], conformal exponent 0, which 
has logarithmic singularities at the punctures, with the 
initially C* field x = exp(— 4(/)). This new variable, along 
with the other BSSN variables, will remain finite pro- 
vided that one uses a suitable choice for the gauge. We 
obtained accurate, convergent waveforms by evolving this 
system in conjunction with a modified 1-1- log lapse, a 
modified Gamma-driver shift condition 0, |43|, and an 
initial lapse set to a = 2x/(l + x)- The lapse and shift 
are evolved with 

{dt-f3'd^)a = -2aK, (5) 
dtP" = B", dtB" = 3/49tf " - r/B°. (6) 

In this paper we will explore how the value of the gauge 
parameter ry affects the quasi-local calculation of the mo- 
mentum. We use the Carpet [4^ mesh refinement 
driver to provide a 'moving boxes' style mesh refinement. 
In this approach refined grids of fixed size are arranged 
about the coordinate centers of both holes. The Car- 
pet code then moves these fine grids about the compu- 
tational domain by following the trajectories of the two 
black holes. We track the location of the app arent hori- 
zons using the AHFinderDirect thorn [i^. 

For this calculation we take in Eq. ([3]) to be (1, 0, 0), 
(0, 1,0), and (0,0, 1). We make this choice regardless of 
slicing. This choice works for the conformally flat Bowen 
and York initial data as noted above. 



TABLE I: Initial data parameters for quasi-circular orbital 
frequency uj/M = 0.05 according to 3PN. In this run we take 
the horizon mass ratio of the holes q = jm^l — 3/8, and 
vanishing spins. Si. The punctures are located along the x- 
axis with momenta P along the y-axis. Puncture mass pa- 
rameters are denoted by mf , and horizon mass by . 





1.7604572 


m\/M 
ml/M 


0.718534207968 


XI 1 


-4.7455652 


0.257487827988 


SUM' 
SUM'' 


0.0000000 


m? /M 


0.735380191 


0.0000000 


/M 


0.27582649 


P/M 
U /M' 


0.10682112 


Madm/M 
J/M' 


1.00001 


0.69498063 


0.69498063 



A. Orbital Linear Momentum 

We consider initially non-spinning, unequal-mass, or- 
biting black-hole binaries, with mass ratio 3 : 8, starting 
from an initial separation that leads to at least two or- 
bits prior to merger. We apply our formula ([3]) to mea- 
sure the linear momentum of each hole as they orbit each 
other. The initial data parameters for this configuration 
are summarized in Table HI We evolved this dataset us- 
ing 10 levels of refinement with the coarsest resolution of 
h = QAM and outer boundary at 320M, and finest res- 
olution oi h = M/80. In addition we evolved the same 
dataset after refining the resolution at each level by a 
factor of 4/5 and 4/6. 

We plot the x and y components of the linear mo- 
mentum of each horizon, as well as the Euclidean norm, 
versus time in Fig. [TJ We have calculated the momen- 
tum both using formula ([3]) and the purely coordinate 
momentum Pi — rriidyi/dt, where iji is the trajectory of 
puncture i and is the horizon mass. Note that the 
coordinate momentum is initially zero due to our choice 
P'it = 0) = B'{t = 0) ^ 0, and that it is consistently 
lower than the initial momentum of the holes (and de- 
creases at late-times). Therefore, we expect that the 
quasi-local momentum provides a more accurate mea- 
surement than the coordinate momentum. Furthermore, 
we expect that the momentum will increase as the binary 
inspirals, in qualitative agreement with the behavior of 
the quasi-local momentum. Note that both evaluations 
of the momentum agree in phase, hence we expect that 
the quasi-local momentum, which has a more accurate 
amplitude, will provide a more accurate measurement of 
the instantaneous angular velocity. 



B. Recoil momentum 

In this section we study the recoil velocity of an already 
merged black-hole binary that acquires linear momen- 
tum as a reaction to the asymmetric radiation of gravita- 
tional waves during the inspiral and merger phases. We 
apply formula ([3]) to a single, slightly distorted, black 
hole and we attempt to compute speeds of the order 
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20 40 60 80 100 120 140 160 

t/M 

FIG. 1: The x and y components and magnitude of the in- 
dividual horizon momentum computed for the unequal mass 
black-hole binary with mass ratio q — 3/8 until merger at 
t ~ 164M. Here P denotes the quasi-local momentum and 
V = dy/dt. 



TABLE II: Initial data parameters for the equal mass opposite 
spinning black holes. The punctures are located along the x- 
axis starting from rest. Spins are denoted by Si, puncture 
mass parameters by , and horizon mass by mf . 



xi/M 


-3.5000 


m?JM 


0.427644 


X2/M 

Sf/M^ 


+3.5000 


0.427644 


+0.15 


/M 


0.517407 


-0.15 


m§ /M 


0.517407 


P/M 
L^/M^ 


0.0000 


Maom/M 
J/M^ 


0.999998 


0.0000 


0.0000 



v/c ^ 0.0001, which is three orders of magnitude smaller 
than the speeds in the orbital case. 

We evolve two equal-mass black holes with opposite 
spins (pointing perpendicular to the line joining the 
holes), starting from rest, at an intermediate separa- 
tion. The initial data parameters are summarized in 
Table |lTl For this configuration the recoil momentum, 
as determined by an extrapolation of ^4 to r = oo, is 
Py/M = (20.4 ± 0.5)km s'^ (P^ = = 0). We evolved 
these data using 9 levels of refinement with coarsest res- 
olution h = 8M and finest resolution h = M/32, as well 
higher resolutions runs with finest resolutions of M/48 
and M/64 (with a corresponding increase in gridpoints 
on all levels). Here we explore the dependence of the 
measured kick both on resolution and the rj parameter 
in Eq. (Note that Eq. ([5]) approximates maximal 

slicing, i.e. X — > 0, at late times.) 

The left panel of Fig. [2] shows the quasi-local recoil 
momentum versus time for rj — I and the three reso- 
lutions. Note the rapid convergence of the asymptotic 
value of the momentum versus resolution. The measured 
convergence rate is greater than 5, but even at a cen- 
tral resolution of M/48 the error in the recoil is within 
3km s~^. An extrapolation to infinite resolution gives 
Py/M = (21.5±0.5)km s~^. The small disagreement be- 
tween the quasi-local recoil and recoil calculated from ■04 




FIG. 2: The remnant recoil velocity for the head-on test as 
determined by the isolated horizon formula. The panel on 
the left shows Py versus resolution for = 1. Note the rapid 
convergence of the asymptotic plateau. The panel on the right 
shows Py versus rj for h = A'f/32. The expected value of the 
recoil is Py/M = (20.4 ± 0.5)km s~\ 




180 190 200 210 220 230 240 

t/M 



FIG. 3: The recoil velocity of the remnant (bottom) and error 
(top) for the orbiting binary versus resolution for rj = 2. 



is due to the use of a non-vanishing rj. The right panel of 
Fig. [5] shows the dependence of the quasi- local momen- 
tum on the gauge via variations in r/. Note that while 
1] has essentially no effect on the recoil calculated from 
ijji, the distortion in the gauge caused by large rj has a 
significant affect on the quasi- local momentum for rj > 2. 
In particular, for large ry, the quasi-local momentum dis- 
plays a slow secular decay towards a final asymptotic 
value, with both an increasing amplitude and decreasing 
decay rate, as rj is increased. 

Note that the quasi-local formula provides an accurate 
measurement of the recoil at i ~ BOM, while a measure- 
ment based on ^'4 requires evolutions of at least 150M 
for observers at 40M. 

Figure [3] shows the quasi-local recoil momentum calcu- 
lation for the orbiting binary configuration as a function 
of time and resolution. The expected recoil velocity in 
this case is y = 175km s~^. 



IV. CONCLUSION 

There is currently no other fully relativistic method to 
compute the linear momentum of black holes in a binary. 
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Knowledge of the linear momentum of each hole provides 
an important diagnostic in comparing fully-nonlinear re- 
sults with post-Newtonian predictions of the trajectory 
and waveform (for instance as an alternative measure of 
the instantaneous angular velocity of the binary system). 

The quasi-local approach we propose represents not 
only an alternative measure of the linear momentum, but 
also provides an accurate measurement of the recoil much 
sooner than can be obtained from the waveform. 

In future work we will study a more coordinate in- 
dependent and robust derivation of isolated horizon in- 
spired formulae to evaluate the linear momentum of black 
holes. This can be achieved by an improved evaluation 
of the in addition to extrapolation to the rj ^ and 
/i — > limits, and making greater use of the horizon ge- 



ometry than was done here. 
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